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FORECASTER,  A  MARKOVIAN  MODEL  TO  ANALYZE 
THE  DISTRIBUTION  OF  NAVAL  OFFICERS 

by 
Paul  R.  Milch 

1.  INTRODUCTION  AND  BACKGROUND 

In  an  earlier  report  [4]  the  author  constructed  an  analytical  model  based 
on  Semi-Markov  Processes  representing  the  career  paths  of  military  officers 
in  general  and  Surface  Warfare  Officers  (SWO's)  of  the  U.S.  Navy  in 
particular.  The  model  was  based  on  a  rectangular  grid  representation  of  the 
entire  career  structure  of  a  group  of  officers.  In  the  rectangular  grid,  rows 
represent  the  various  activity  types  (e.g.,  sea  duties  and  shore  duties)  in 
which  an  officer  of  the  group  whose  career  structure  is  modeled  may  be 
engaged  from  time  to  time.  Columns  of  the  grid  stand  for  the  succession  of 
tours  that  makes  up  the  career  of  any  military  officer.  In  this  way,  an  officer's 
career  may  be  represented  by  a  path  through  a  number  of  nodes  (billets) 
whose  first  coordinates  show  the  activity  types  and  the  second  coordinates  the 
tour  numbers. 

Associated  with  each  node  on  the  grid  is  a  specific  non-negative  integer, 
called  the  tourlength,  which  determines  the  length  of  time  the  officer  spends 
in  that  node  or  billet.  A  zero  tourlength  means  that  the  node  or  billet  is 
infeasible.  This  way,  a  matrix  of  size  A  by  R  determines  the  tourlength  matrix 
where  A  and  R  are,  respectively,  the  number  of  activity  types  and  the  number 
of  tours  in  the  entire  career  structure.  A  similar  A  by  R  matrix  may  represent 
the  number  of  officers,  called  incumbents,  currently  occupying  billets  in  all 


the  feasible  nodes.  Yet  a  third  A  by  R  matrix  of  non-negative  integers 
provides  the  numbers  of  accessions  or  recruits  entering  the  system  at  future 
times.  Typically,  this  latter  matrix  is  all  zeros,  except  for  some  elements  of  the 
first  column,  because  new  entries  into  the  system  should  start  with  the  first 
tour.  However,  in  some  officer  communities,  such  as  the  Navy  Nurse  corps, 
for  example,  there  are  significant  numbers  of  "recalls"  and  other  "lateral 
entries"  who  join  the  community  in  later  tours.  Finally,  the  progression  of 
advancement  on  the  career  grid  is  governed  by  a  sequence  of  transition 
probability  matrices  which  prescribe  the  probability  of  an  officer  currently  in  a 
specific  node  (i.e.,  activity-tour  combination)  transiting  to  another  (possibly 
the  same)  activity  type  upon  finishing  his/her  current  tour.  These  matrices 
are  all  A  by  A  and  there  are  one  fewer  of  them  than  the  total  number  of  tours 
(i.e.,  R-l).   A  more  detailed  description  of  the  system  is  given  in  Milch  [4]. 

In  the  same  report  the  mathematical  details  and  formulas  are  also 
described  for  the  expected  numbers  of  officers  occupying  each  node  (billet)  on 
the  grid  at  some  future  point  in  time.  This  is  based  on  first  computing 
probabilities  that  an  officer  occupying  a  node  (billet)  now  will  be  anywhere  on 
the  grid  t  time  units  later.  This  computation  is  based  on  a  Semi-Markov 
Process  representation  of  the  position  occupied  by  an  officer  at  time  t.  When 
the  appropriate  formulas  were  programmed  for  computation  in  APL  on  the 
IBM  3033  mainframe  computer  at  the  NPS  it  was  found  that  the  computation 
required  either  excessively  large  amounts  of  memory  or  huge  amounts  of 
CPU  time.  This  occurred  even  after  various  time-saving  measures  were 
effected.  Since  it  was  desirable  to  make  the  model  available  for  personal 
computers  as  well,  the  necessity  of  revising  the  computational  core  became 
even  more  acute. 


Simultaneously  with  efforts  to  revise  the  mathematical  core  of  the  model 
two  attempts  have  been  made  to  evaluate  the  practical  usefulness  of  the 
model  apart  from  its  computational  speed  (or  rather  lack  of  it).  Thus  Johnson 
[3]  built  an  interactive  user  friendly  interface  for  the  model  now  labeled 
FORECASTER,  thereby  greatly  enhancing  its  usability.  He  also  demonstrated 
how  FORECASTER  may  be  used  by  analyzing  the  currently  topical  issue  of 
joint  duty  assignments  effected  by  the  recently  passed  Goldwater-Nichols 
Department  of  Defense  Reorganization  Act.  Title  IV  of  this  act  mandates  the 
creation  of  so-called  joint  duty  officers  in  all  four  military  services  and 
prescribes  their  professional  education,  as  well  as  their  use  by  the  services. 
Johnson  examined  the  effect  of  the  new  law  on  the  SWO  community:  its 
ability  to  fill  its  "fair  share"  of  joint  duty  billets  and  how  that  may  affect  the 
traditional  SWO  career  path.  More  recently,  Drescher  [11  examined  the  same 
problem  in  another  Navy  officer  community,  namely  that  of  Tactical 
Aviation  pilots  and  flight  officers.  At  the  same  time,  Drescher  reorganized 
the  APL  functions  used  in  FORECASTER  and  provided  a  more  thorough 
documentation  for  the  model. 

Concurrently  with  Drescher's  work,  this  author  has  revised  the 
mathematical  core  of  the  model.  This  new  mathematical  formulation,  a 
conceptually  much  simpler  one  than  the  original  presented  in  Milch  [4] 
turned  out  to  be  significantly  faster  computationally  as  well  and  that  is  the 
subject  of  this  report.  It  is  also  hoped  that  the  new  mathematical  formulation 
of  the  model  will  make  it  easier  to  effect  another  significant  improvement  in 
the  usefulness  of  the  model.  This  issue  has  to  do  with  the  way  the  user  must 
input  data  into  the  current  model,  as  well  as  the  way  the  user  must  make  use 
of  the  output  results. 


Briefly  summarized,  this  aspect  is  due  to  the  fact  that  the  model  is  based 
on  a  sequence  of  tours  that  the  officers  undergo.  This  is  a  realistic  simulation 
of  what  actually  happens,  but  it  also  means  that  incumbents  data  (i.e., 
numbers  of  officers  currently  occupying  billets)  must  be  categorized  by 
activities  and  tours.  Unfortunately,  the  usual  method  of  recordkeeping  of 
officers  in  any  Navy  community  is  not  organized  by  tours,  but  rather  by  grade 
and  years  of  commissioned  service.  Likewise,  results  of  the  model  categorize 
officers  distributed  at  some  future  time  by  activity  and  tour  number.  These 
results  must  then  be  compared  to  billet  requirements  that  are  usually  not 
organized  by  tour  number  but  by  grade. 

Johnson  [3]  and  Drescher  [1]  demonstrated  that  the  data  is  available  by 
tour  number  from  the  source  where  all  officer  data  originates,  the  Officer 
Master  File.  However,  the  fact  remains  that  standard  personnel  planning  by 
all  Navy  offices  where  such  work  is  done  is  based  on  officer  and  billet  data 
that  are  not  organized  by  tour  numbers  and  this  practice  is  not  likely  to 
change  in  the  near  future.  For  this  reason  it  is  desirable  to  convert  the  model 
to  accept  input  and  produce  output  by  grade  and /or  years  of  service  rather 
than  tour  number  in  order  to  accommodate  user  needs. 

This  effort  is  currently  under  way.  In  the  meantime,  the  new 
mathematical  core  was  designed  to  accept  data  as  the  original  model  did 
which  is  described  in  Milch  [4]  and  outlined  at  the  beginning  of  this  section. 
The  next  section  describes  the  mathematical  details  of  the  model  according  to 
the  newly  created  APL  code  available  on  either  the  IBM  3033  mainframe 
computer  or  a  personal  computer. 


2.     DATA  PREPARATION 

As  a  point  of  departure,  the  notation  of  the  author's  previous  report  [4] 
will  be  adopted  wherever  appropriate.  Thus,  the  tourlength  matrix,  L,  has  the 
elements 

duration  of  time  spent  in  activity  type  i 
*in  ~  by  an  officer  on  his/her  nth  tour, 

for    i  =  1, ...,  A    and    n  =  1, ...,  R. 

The  incumbents  matrix,  Q,  is  defined  as  having  elements 

number  of  officers  at  time  zero  occupying 
^in~  activity  type  i  billets  during  their  nth  tour, 

for    i  =  1, ...,  A    and     n  =  1, ...,  R. 

Further,  the  accessions  matrix,  C,  is  defined  as  consisting  of  elements 

number  of  officers  accessed  during  a  time  period 
Cin  ~  directly  into  an  activity  type  i  and  tour  number  n, 

for    i  =  1, ...,  A    and    n  =  1, ...,  R. 

Most  of  the  time,  matrix  C  will  have  positive  elements  only  in  its  first 
column  (for  n  =  1),  since  a  newly  accessed  officer  must  usually  start  in  tour 
number  one.  It  is  interesting  to  note,  however,  that  this  is  not  always  the 
case.  For  example,  officers  are  sometimes  recalled  shortly  after  their 
separation  from  the  service  or  are  laterally  transferred  from  one  community 
to  another  and  may  thus  join  the  community  being  analyzed  in  a  tour  other 
than  their  first  one. 

Finally,  the  transition  probability  matrices,  P(n)  are  defined  to  have 
elements, 


probability  that  an  officer  transfers  to  an  activity  type  j 
pij(n)    =    on  his/her  n+lst  tour  upon  completion  of  his/her 
nth  tour  in  an  activity  type  i, 

for    i,j  =  1, ...,  A    and    n  =  1, ...,  R-l. 

The  construction  of  the  model  will  be  described  in  steps  like  those  of  an 

algorithm  as  that  is  most  suitable  for  the  procedure  which  turned  out  to  be 

both  simple  and  computationally  efficient. 

Step  1:      Unravel  the  tourlength  matrix,  L,  by  columns,  starting  with  the 

A 

first  column,  followed  by  the  second,  etc.,  into  a  vector  £  .    More 
precisely,  define 

A         A  A 

*=(*i,...,/k) 

where 

A 

^k  =  4n    when    k  =  (n-l)A+i 
for    i  =  1, ...,  A    and    n  =  1, ...,  R-l    so  that    k  =  1,  2, ...,  K  =  AR 
Step  2:      Define  the  indices  k],  ki, ...,  kj<*  as  the  first,  second,  ...,  last  values 

A 

of  k  for  which  i^  is  positive.    These  are  called  the  "ascending 
ladder  indices  or  epochs"  by  Feller    [2].    Then  the  abbreviated 

A 

version  of  £  may  be  defined  by  omitting  its  zero  components.    In 
other  words, 


t*  =  Ukl *kK.) 


where 


R      A 

K*=n  Sin    With    8in= 
n=l    i=l 


1     if    4n>0 
0     if    4n=0 


It  will  be  useful  to  define  the  partial  sums  of  the  components  of 


1  K* 


as  well.  Let 

u  u     A 

so  =  0    and    su  =  X  **  =X  t  k  f°r  u  =  1, ...,  K* 
a  j=i    J     j=i      ' 

since      /*  =  £  k .    In  Feller's  terminology  the  su  are  the  "ascending  ladder 

heights." 

Step  3:  Next  the  vector  /*  is  expanded  by  repeating  each  of  its 
components  exactly  as  many  times  as  indicated  by  its  value.  That 
is,  if  £*  =  (3,  1,  5)  then  the  expanded  vector  £**  will  have  nine 
components  with  the  first  three  equal  to  3,  the  fourth  equal  to  1 
and  the  last  five  equal  to  5.  In  general, 

1  K** 


where 


R       A 
n=l     i=l 


The  components  of  £**  may  formally  be  defined  as 


£     =  £       when    su_i  +  l<j<  su  and    u  =  1, ...,  K*. 

J  ku 

For  example,  when  u  =  1,  /     =  £     for  all  j  =  1, ...,  s\  =£    .    This  means  that 

j        ki  ki 

for  the  above  example  /     =  3,  for  j  =  1,  2,  3. 

j 

Step  4:      Next   a   similar   procedure   is   followed    with   respect   to    the 
incumbents  matrix,  Q.    First,  Q  is  unraveled  by  columns  into  the 

A 

vector  q.  That  is, 


A  A  A 


g  =  (qr-vqk) 

where 

qk  =  qin    when    k  =  (n-l)A+i 

for    i  =  1, ...,  A    and    n  =  1, ...,  R    so  that    k  =  1, ...,  K  =  AR. 

Then  cj.  is  abbreviated  by  the  omission  of  those  of  its  components,  qk,  for 

A 

which  the  corresponding  I  k  is  zero.  That  is, 


a*  =  (q*/  •••/  q*  ) 

1  K 


*        A 


where    q  =  q,     for   j  =  1, ...,  K*. 
J         ' 
Note  that  the  decision  to  keep  or  omit  a  component,  qk,  of  q  is  based  not 

on  whether  qk>0  or  qk=0,  but  on  whether  /k>0  or  ^=0- 

Finally,  the  vector  cj*  is  expanded  to  repeat  each  of  its  components  as 

many  times  as  indicated  by  the  value  of  the  corresponding  component  of  /*. 

More  formally, 


jT-Cq^-vq      > 


where 


**■ 


q     =q       when    su_i  -t- 1  <j<su  and    u  =  1, ...,  K* 

j         K 


Step  5:      Next  two  "flow  rate"  vectors  are  defined  by  dividing  the  number 
of  incumbents  by  their  respective  tourlengths: 

r*  =  (r*, ...,  r*  )    and    r**  =  (r**,  ...,  r**  ) 


8 


where 


r    =— *■ 

j  e 


for    j  =  l,...,K* 


and 


** 

q. 

r**=  -1-    for    j  =  l, ...,  K' 


Step  6:       A  slightly  different  procedure  is  followed  with  the  accessions 
matrix,  C.   Again,  C  is  first  unraveled  by  columns  into  the  vector 


where 


A  A  A 

c  =(ci,  ...,ck) 


Ck  =  Qn    when    k  =  (n-l)A+i 
for    i  =  1, ...,  A    and    n  =  1, ...,  R    so  that    k  =  1, ...,  K  =  AR. 

A 

Then  c  is  abbreviated  by  the  omission  of  those  of  its  components,  c  ^  for 

A 

which  the  corresponding  (.^  is  zero.  That  is, 


A 


where 


c*  =  (c  , ...,  c    ) 

1  K* 


=  ck)    for   j  =  l, ...,  K* 


Then  c*  is  expanded  to  have  K**  components,  but  instead  of  repeating  its 
components  (as  was  the  case  in  expanding  q*  into  cj**)  the  extra  components 
in  c**  are  simply  zeroes.   More  precisely, 


c**  =  (c    ,  ...,  c      ) 

1  K** 


where 


c 


* 
*»_Jc         when    j  =  su-i  +  1    and    u  =  1, ...,  K* 


k 

u 


^0         otherwise. 


These  six  steps  accomplish  the  task  of  rearranging  the  data  originally 
organized  for  user  convenience  into  formats  that  fit  the  need  of  efficient 
storage  and  computation.  Namely,  the  tourlength  matrix,  L,  was 
manipulated  first  to  eliminate  infeasible  billet  types  and  then  to  repeat  its 
positive  components  in  preparation  of  the  computation  of  the  flow  rate 
vectors  r*  and  r**.  Similarly,  the  incumbents  matrix  was  first  transformed 
into  a  vector  in  which  any  incumbents  occupying  infeasible  billets  (due  to  a 
data  recording  error  or  for  any  other  reasons)  have  been  omitted.  The 
practical  effect  of  this  is  that  such  personnel  will  be  attrited  by  the  model 
during  the  forecasting  process  which  is  probably  the  appropriate 
representation  of  what  happens  to  such  personnel,  anyway. 

The  purpose  of  dividing  the  number  of  incumbents  by  the  corresponding 
tourlength  and  then  repeating  that  quotient  as  many  times  as  indicated  by  the 
value  of  the  tourlength  is  to  distribute  personnel  occupying  a  billet  type 
evenly  in  terms  of  longevity  (or  experience)  in  that  billet.  This  may  not  be 
the  same  as  the  actual  distribution  in  reality,  but  the  latter  may  either  not  be 
available  or  would  present  too  much  trouble  for  the  user  to  enter  into  the 
model  for  the  additional  accuracy  that  would  be  gained  thereby. 

The  accessions  matrix  is  similarly  treated  except  that  instead  of  repeating 
its  components  all  accessions  are  placed  into  the  first  (i.e.,  zero)  longevity  (or 
experience)  category  and  the  additional  components  of  c**  are  simply  zeroes. 


10 


The  next  task  is  to  consolidate  the  transition  probability  matrices,  P(n),  n  = 
1,  ...,  R-l,  into  a  single  large  transition  probability  matrix  that  will  affect  the 
transitions  of  personnel  occurring  each  unit  of  time.  Conceptually  the 
simplest  way  of  accomplishing  this  would  be  to  construct  a  Markov  Chain 
with  K**  states  corresponding  to  the  components  of  the  vector  r**.  Then 
each  unit  of  time  components  of  r**  would  be  transferred  in  accordance  with 
a  K**  by  K**  transition  probability  matrix  which  would  contain  all  the 
positive  elements  of  all  the  P(n)  matrices,  but  would  otherwise  consist  of 
mostly  zeroes  and  many  ones,  representing  the  fact  that  personnel  cannot 
move  backward  (or  "sideways")  in  time  and  once  placed  in  a  billet  personnel 
must,  by  necessity,  move  up  by  one  longevity  (experience)  level  per  unit  of 
time  until  their  experience  level  equals  the  tourlength  in  that  billet  at  which 
point  the  appropriate  elements  of  the  relevant  P(n)  matrix  would  distribute 
them  into  various  activities  of  their  next  tour. 

The  problem  with  this  approach  is  that  K**,  being  the  sum  of  all 
tourlengths,  is  in  the  range  of  about  300  or  more  for  most  problems  for  which 
this  model  is  being  designed  when  tourlengths  are  given  in  terms  of  quarters 
(3  months)  which  is  probably  the  preferred  unit  of  time,  with  the  alternative 
being  months  (which  in  turn  would  further  increase  the  number  of  states 
three-fold).  When  forecasting  several  units  of  time  into  the  future  having 
the  number  of  states  in  the  neighborhood  of  300  would  tax  the  computational 
capabilities  of  most  micro-computers  in  terms  of  memory  and/or  CPU  time 
and  perhaps  even  of  most  commonly  available  mainframe  computers  for 
which  this  model  could  be  designed.  Therefore,  practicality  dictates  a  better 
approach  which  is  less  consumptive  of  CPU  time  as  well  as  memory. 
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The  problem  of  too  high  dimensionality  will  be  solved  by  using  the 
Markov  Chain  approach  not  to  forecast  the  vector  r**  but  the  much  lower 
dimensional  vector  r*.  This  reduces  the  dimensions  from  K** 
(approximately  300)  to  K*  (approximately  50  or  60).  This  requires  the 
transition  probability  matrix  to  be  of  only  K*  by  K*. 

It  is  this  matrix  that  is  constructed  next  from  the  sequence  of  transition 
probability  matrices  P(n),  n  =  1, ...,  R-l.  This  is  accomplished  in  several  steps. 

First  the  matrices,  P(n),  are  made  consistent  with  the  information 
contained  in  the  tourlength  matrix,  L.  Namely,  if  the  element  lin  is  zero, 
indicating  that  activity  type  i  for  tour  number  n  is  not  feasible  then  the  ith  row 
of  the  matrix  P(n)  is  superfluous  and  in  fact  should  contain  all  zeroes.  If 
some  elements,  PjXn),  are  positive  for  some  j  =  1, ...,  A  it  must  be  due  to  some 
error,  most  likely  an  omission  to  change  these  elements  to  zero  by  a  user  of 
the  model  when  a  formerly  feasible  billet  was  made  infeasible.  At  any  rate, 
/in  =  0  implies  that  Pjj(n)  should  be  zero  for  all  j  =  1,  ...,  A  and  therefore  these 
elements  occupy  redundant  space  in  P(n).  Similarly,  the  ith  column  of  the 
matrix  P(n-l)  should  contain  all  zero  elements  and  is  likewise  redundant. 
This  results  in  the  following 

A 

Step  1:      Define  P(n),  for  n  =  1,  ...,  R-l,  to  have  elements 


P„(n)  = 


p..(n)      if     I.  I,    ,  >0 
r  lj  in  jn+1 


for  i,j  =  1,  ...,  A. 


10  if      Wjn+1    =0 

A 

Then  define  the  matrix  P*(n),  for  n  =  1,  ...,  R-l,  to  be  the  same  as  P(n)  except 
that  any  all-zero  row  or  all-zero  column  is  omitted.  Note  that  the  matrices, 
P*(n),  may  no  longer  be  square  matrices  and  in  fact  it  will  be  useful  to 
introduce  notation  for  their  row  and  column  dimensions.    Namely,  let 
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dr(n)  =  no.  of  rows  of  P*(n) 
dc(n)  =  no.  of  columns  of  P*(n) 

for  n  =  1, ...,  R-l. 

Also,  let  dr(0)  =  dc(0)  =  0. 

Further,  define  the  partial  sums  of  these  row  and  column  dimensions: 

n 

Dr(n)  =  X  dr(m) 

m=0 
n 

Dc(n)  =  X  dc(m> 

for  n  =  0,  1, ...,  R-l. 
Note  that 

A  A 

dr(n)  =  X  ^n    and    dc(n)  =  X  Sin+i 
i=l  i=l 

forn  =  1, ...,  R-l. 

Therefore,  dc(n)  =  dr(n+l)  or  the  column  dimension  of  P*(n)  is  the  same 
as  the  row  dimension  of  P*(n+1). 

Also,  note  that 

R     A  R-l  A 

K*  =  £  X  Sin  =  X  dr(n)  +  X  &iR  =  Dr(R-l)+dc(R-l) 
n=l   i=l  n=l  i=l 

and  similarly, 

R-l   A  A  R-i 

K*=n  5in+i  =X  8ii  +  X  dc(n)  =  dr(l)+Dc(R-l) 
n=0  i=l  i=l  n=l 

Therefore, 

Dr(R-l)  +  dc(R-l)  =  dr(l)  +  Dc(R-l) 
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or 

Dr(R-l)  -  dr(l)  =  Dc(R-l)  -  dc(R-l) 

That  is,  the  sum  of  all  but  the  first  row  dimensions  equals  the  sum  of  all  but 
the  last  column  dimensions. 

Step  2:      The  K*  by  K*  dimensional  transition  probability  matrix,  P*,  may 
now  be  introduced  with  all  elements  zero  except  the  following: 


p  =  p  (n) 

rDr(n-l)+i,dr(l)+Dc(n-l)+j     Mj 


for    i  =  1, ...,  dr(n)    and    j  =  1, ...,  dc(n)    and    n  =  1, ...,  R-l. 

Note  that  by  necessity  the  first  dr(l)  columns  of  the  matrix,  P*,  must  be  all 
zeroes.  Similarly,  the  last  dc(R-l)  rows  of  the  matrix  P*,  are  also  all  zeroes.  It 
is  also  clear  that  P*  is  an  upper  triangular  matrix,  in  fact  all  its  main  diagonal 
elements  are  also  zero. 

The  following  example  may  serve  as  an  illustration  of  the  entire 
procedure  described  above.  The  three  matrices,  L,  Q  and  C  are  given  as 


L  = 


( 


V 


and 


with    i  =  3    and    R  =  4. 


2 

0 

3 

0 

\ 

( 

50 

0 

60 

0 

A 

0 

1 

0 

2 

Q  = 

0 

10 

8 

50 

3 

2 

0 

4 

) 

\ 

90 

30 

20 

80 

J 

( 

7( 

)      0 

10 

0 

A 

C  = 

\ 

1( 

5( 

)      0 
)       0 

5 
0 

0 
0 

) 

Then  the  unraveling  process  produces  the   K=12  dimensional  vectors 


I  =(2,0,3,0,1,2,3,0,0,0,2,4) 
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g  =  (50,  0,  90, 0,10,  30,  60,  8,  20,  0,  50,  80) 


c  =  (70,10,  50, 0, 0, 0, 10, 5, 0, 0, 0, 0) 

A 

Examining  t  shows  that  the  "ascending  ladder  epochs"  are: 

ki  =  l,  k2  =  3,   k3  =  5,   k4  =  6,   k5  =  7,   1^  =  11,   k7  =  12. 
Therefore,  the  abbreviated  vectors  are  K*  =  7  dimensional,  namely: 

/*=  (2,  3, 1,2, 3, 2, 4) 
g*  =  (50, 90, 10, 30,  60, 50,  80) 
c*  =  (70,  50, 0, 0, 10, 0,  0) 
The  "ascending  ladder  heights"  are: 

So  =  0,    Si=2,    S2  =  5,    S3  =  6,    54  =  8,55  =  11,    S6  =  13,     S7  =  17. 

The  vectors  with  repeated  components  are  K**  =  17  dimensional,  namely: 

C*  =  (2,  2,  3,  3,  3, 1,  2,  2,  3,  3,  3,  2,  2,  4,  4,  4,  4) 

a,**  =  (50,  50, 90, 90, 90, 10, 30, 30,  60,  60,  60, 50, 50,  80, 80,  80,  80) 

c**  =  (70,  0,  50,  0,  0,  0, 0, 0, 10, 0,  0, 0,  0, 0,  0,  0,  0) 
Next,  the  two  flow  rate  vectors  of  dimensions  K*  =  7  and  K**  =  17  are 
computed  as: 

r*  =  (25,  30, 10, 15,  20,  25,  20) 

r**  =  (25,  25,  30,  30,  30, 10, 15, 15,  20,  20, 20,  25,  25,  20,  20,  20,  20) 
Assume  that  the  following  three  P(n),  n  =  1,  2,  3,  matrices  are  also  given: 
(  0     .5     .4  A  (  .8     0     .1  *\  (  .2    .6     .1  ^ 


P(l)  = 


.1     .7     0 
,•0     -8     -1 


P(2)  = 


.8     0      0 
^.7     .1     0 


P(3)  = 


.2     .5     .2 
^.8     .1     0 
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The  first  step  is  to  replace  any  positive  pij(n)'s  with  zeroes  if  either  £[n  =  0  or 


/jn=0. 


A 
P(D  = 


(  0     .5     A\ 

0     0     0 
,0    .8     .1 


P(2)  = 


(  o    o   o  A 

.8     0     0 
,•700, 


P(D  = 


(  0     .6     .1   ^ 

0     0      0 
^  0     0      0    J 


Next,  all-zero  rows  and  all-zero  columns  are  eliminated: 


P*(D  = 


(  -5     .4  ^ 


.8     .1 


P*(2)  = 


f.8\ 


P*(3)  =  (.6    .1) 


Note  that 

dr(l)  =  2,  dr(2)  =  2,  dr(3)  =  1  and  dc(l)  =  2,  dc(2)  =  1,  dc(3)  =  2. 
Therefore, 

Dr(0)  =  0,  Dr(l)  =  2,  Dr(2)  =  4,  Dr(3)  =  5 

and  Dc(0)  =  0,  Dc(l)  =  2,  Dc(2)  =  3,  Dc(3)  =  5. 

The  K*=7  dimensional  square  matrix  P*  is 


P*  = 


(  ° 

0 

(.5 

V.8 

:?) 

0 

0 

0  \ 

0 

0 

0 

0 

0 

0 
0 

0 
0 

0 
0 

0 
0 

(:?) 

0 
0 

0 
0 

0 

0 

0 

0 

0 

(.6 

.1) 

0 

0 

0 

0 

0 

0 

0 

\° 

0 

0 

0 

0 

0 

°) 

Note  the  way  the  matrices,  P*(l),  P*(2)  and  P*(3),  are  embedded  in  the 
large  P*  matrix,  in  the  form  of  a  "descending  staircase"  which  is  typical  of  the 
structure  of  the  P*  matrix. 

3.     FORECASTING 

Forecasting  will  be  accomplished  in  several  steps  using  both  K*-  and  K**- 
dimensional  vectors. 
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The  first  five  of  these  steps  are  taken  for  a  single  value  of  t,  starting  with 
t  =  0  and  then  repeated  for  t  =  1, ...,  T-l  where  T  is  the  forecasting  horizon 
for  which  the  distribution  of  officers  is  required.  After  the  first  five  steps 
have  been  successfully  repeated  for  t  =  0,  1,  ...,  T-l  the  actual  forecasting 
procedure  is  complete.  Steps  6,  7,  and  8  then  serve  the  purpose  of  rearranging 
the  result  in  the  standard  A  by  R  matrix  format. 

Step  1:  First  the  K*-dimensional  vector  r*(t)  is  established,  made  up  of 
those  elements  of  r**(t)  that  stand  for  the  highest  "experience" 
level  in  each  billet.  More  precisely,  using  the  "ascending  ladder 
height/' 


r*(t)  =  (r(t),...  ,r    (0) 

1  K* 


where 


r*(t)  =  r**(t)  for  u  =  1,  2,  ...,K* 


with    r**(0)  =  r**   as  given  in  the  previous  section. 
Step  2:      Update  r*(t)  by  the  usual  Markovian  equation  of 

r**(t+l)  =  r*(t)P* 

Step  3:      Update    r**(t)    by  "shifting"  its  components  one  component  to 
the  right,  that  is, 


r.    (t+1) 

J 


0  if  j  =  1 


** ,  > 


*  * 


r     (t)  ifj  =  2,...,K** 


with  r**(0)  =  r**  as  in  Step  1. 
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Step  4:  Replace  those  components  of  the  updated  vector,  r**(t+l),  which 
stand  for  zero  "experience"  level  in  each  billet  by  successive 
components  of  the  updated  vector  r*(t+l).   More  precisely, 

r**     (t+l)  =  r*     (t+1)    for    u  =  0,  ...,KM. 
su+l  u+1 


Step  5:  As  a  last  step  in  the  updating  process  the  accession  vector,  c**  is 
added  to  the  vector  r_**(t+l)  obtained  in  Step  4.  That  is, 

new  r**(t+l)  =  old  r**(t+l)  +  c  ** 

These  five  steps  are  now  repeated  for  the  next  higher  t-value  until  r**(T) 
is  obtained.  Then  Steps  6,  7,  and  8  follow: 

Step  6:  To  add  the  appropriate  components  of  r  **(T)  is  easiest 
accomplished  in  terms  of  the  "ascending  ladder  heights,"  su. 
Namely,  define 

SjfCI)  =  (cfCT), ...,  q*Cn)    for    T  =  l,2,... 


where 


q*    JT)=    I  r"(T)    for    u  =  0, 1, ...,  K*-l 
u+1  j-^+l  ) 


Step  7:  The  vector  cj*(T)  is  expanded  to  insert  zeroes  at  the  appropriate 
places,  representing  infeasible  billets.  This  is  accomplished  in 
terms  of  the  "ascending  ladder  epochs,"  kj.  Define 


g(D  =  (qi(D/...,qK(D) 


where 
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q>=< 


q  (T)  whenk  =  kj    for     j  =  l,...,K* 

0  otherwise. 


Step  8:       Finally,  the  a,  (T)  is  re-raveled  into  the  original  data  format  of  an 
A  by  R  matrix.   That  is,  the  matrix  Q(T)  is  defined  with  elements 

where    k  =  (n-l)A+i,    with    i  =  1, ...,  A    and    n  =  1, ...,  R. 
The  numerical  example  introduced  at  the  end  of  the  previous  section 
is  continued  to  illustrate  the  procedure  outlined  above.     The  forecasting 
horizon  of  T=2  will  be  used. 

In  Step  1,  the  "ascending  ladder  heights" 

so  =  0,  si  =  2,  S2  =  5,  S3  =  6,  S4  =  8,  S5  =  11,  S6  =  13,  S7  =  17 

are  used  to  find  the  K*  =  7  dimensional  vector 

r*(0)  =  (25,  30, 10, 15,  20,  25,  20). 

In  Step  2,  Markovian  forecasting  is  applied  to  this  vector  using  the  7  by  7 
P*-matrix  given  at  the  very  end  of  the  previous  section.   Thus, 
r*(l)  =  r*(0)P*  =  (0,  0,  36.5, 13, 18.5, 12,  2). 

In  Step  3,  the  K**  =  17  dimensional  vector,  r**(0)  =  r**  is  shifted  one 
component  to  the  right  resulting  in 

r**(i)  =  (0,  25,  25,  30,  30,  30, 10, 15, 15,  20,  20,  20,  25,  25,  20,  20,  20). 

Next,  in  Step  4,  using  the  "ascending  ladder  heights,"  again,  appropriate 
components  of  r**(l)  above  are  replaced  by  successive  components  of  r**(l) 
computed  in  Step  2.   The  result  is 
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r**(l)  =  (0,  25,  0,  30,  30,  36.5, 13, 15, 18.5,  20,  20, 12,  25,  2,  20,  20,  20) 

where  the  replaced  components  are  marked  by  an  arrow  above  them. 
In  Step  5,  the  accession  vector 

c**  =  (70,  0,  50, 0,  0,  0,  0,  0, 10,  0,  0,  0,  0,  0,  0,  0,  0) 

is  added  to  the  result  of  Step  4  to  obtain 

r**(l)  =  (70, 25,  50, 30,  30, 36.5, 13, 15,  28.5, 20,  20, 12,  25, 2,  20,  20, 20). 

Next,  Steps  1  through  5  are  repated  for  t  =  1  with  the  following  results: 
Step  1:    r*(l)  =  (25, 30,  36.5, 15,  20,  25, 20) 
Step  2:    r(2)  =  r*(l)P*  =  (0, 0, 36.5, 13, 39.7, 12, 2) 
Step  3:    r**(2)  =  (0,  70, 25,  50,  30, 30, 36.5, 13, 15,  28.5, 20,  20, 12,  25, 2,  20, 20) 

Step  4:     r**(2)  =  (0, 70,  0, 50, 30, 36.5, 13, 13, 39.7, 28.5,  20, 12, 12,  2,  2,  20, 20) 

Step  5:  new  r  **(2)  =  (70,  70,  50,  50,  30, 36.5, 13, 13,  49.7, 28.5,  20, 12, 12,  2,  2,  20,  20) 
Since  r**(T)  has  been  obtained,  it  remains  only  to  reformat  the  result. 
In  Step  6,  the  "ascending  ladder  heights,"  are  used  again  to  compute 

components  of  the  vector 

0/(2)  =  (140,  130,  36.5,  26,  98.2,  24,  44). 

The  details  are  as  follows: 

2  5 

q*(2)=  Xr**(2)  =  70  +  70  =  140         q*(2)=  £r**(2)=  50  +  50  +  30  =  130 
t    j=l  J  ^    j=3  J 

6  8 

q*(2)=Xr**(2)  =  36.5  q*(2)  =  Xr**(2)  =  13  +  13  =  26 

^3    j=6  j  M4    pj  j 

11  13 

Q*(2)=  £r**(2)  =  49.7+28.5  +  20  =  98.2    q*(2)  =  I  r**(2)  =  12  +  12  =  24 
5    j=9  i  *    j=12 ) 


70 


17 

and  q*(2)  =  £  r**(2)  =  2  +  2  +  20  +  20  =  44. 

7  j=14   ) 

In  Step  7,  the  "ascending  ladder  epochs," 

kT  =  1,  k2  =  3,  k3  =  5,  lq  =  6,  k5  =  7,  k*,  =  1 

are  used  to  expand  the  K*  =  7  dimensional  vector  cj.*(2)   into   the  K=12 
dimensional  vector 

oU)  =  (140,  0,  130,  0,  36.5,  26,  98.2,  0,  0,  0,  24,  44). 
Finally,  in  Step  8,  this  vector  is  re-raveled  into  the  3  by  4  matrix 

f  140       0       98.2     0    A 
Q(2)=        0       36.5       0       24 

v  130      26         0       44  J 

In  practice,  each  element  of  this  matrix  would  likely  be  rounded  to  the 
nearest  integer  if  this  is  indeed  the  end  of  the  forecasting  process,  as  these  are 
the  expected  numbers  of  officers  two  quarters  into  the  future. 

4.     AN  APPLICATION. 

A  typical  application  of  the  model  to  the  Surface  Warfare  Officer  (SWO) 
community  of  the  U.S.  Navy  follows.  The  choice  of  activity  types,  listed  in 
Table  1,  was  dictated  by  both  the  nature  of  the  community  and  the  topical 
issue  of  "jointness"  referred  to  in  the  Introduction  and  Background.  Thus,  to 
analyze  the  influence  of  the  new  legislation  on  SWO  career  paths  the  activity 
types  of  "JPME"  (Joint  Professional  Military  Education)  and  "Joint  Tours," 
were  broken  out  separately  in  addition  to  "Postgraduate  Education"  and 
"SWO  Education"  (such  as  Department  Head  School)  from  the  other  "Shore 
Duty,"  with  "Fleet  Unit"  being  the  main  activity  why  the  community  exists. 
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The  tourlength  matrix,  L,  the  Incumbents  matrix,  Q,  and  the  Accessions 
matrix,  C,  are  given  in  Tables  2,  3,  and  4,  respectively.  Finally,  the  eleven 
transition  probability  matrices  are  presented  in  Table  5. 

The  model  was  then  used  to  forecast  the  distribution  of  SWOs  four,  ten 
and  forty  quarters  into  the  future  with  Tables  6,  7  and  8  showing  the  results. 
The  CPU-times  required  to  perform  these  calculations  were,  respectively,  7.69, 
16.31,  and  59.1  seconds  on  a  Zenith  Z-248(R)  microcomputer  with  1.1  MB 
RAM. 

The  model's  user-friendly  features  which  enable  the  user  to  effect  changes 
in  any  of  the  data  input  values  given  in  Tables  1  through  5  and  to  carry  out  a 
billets  versus  officers  comparison  in  a  reasonably  efficient  manner  are 
detailed  in  Johnson  [3]  and  are  not  discussed  here. 
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TABLE  1.  CURRENT  LIST  OF  ACTIVITY  NAMES 
IN  THE  SURFACE  COMMUNITY 


1. 

POSTGR  EDUC 

2. 

JPME 

3. 

JOINT  TOUR 

4. 

SWO  EDUC 

5. 

FLEET  UNIT 

6. 

SHORE  DUTY 

TABLE  2.  CURRENT  TOUR  LENGTHS  IN  QUARTERS 


ACTIVITIES/TOURS 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

1.     POSTGR  EDUC 

0 

8 

8 

8 

6 

6 

6 

6 

6 

0 

6 

0 

2.     JPME 

0 

0 

2 

0 

2 

2 

2 

2 

2 

2 

2 

2 

3.     JOPNT  TOUR 

0 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

4.     SWO  EDUC 

0 

2 

2 

2 

2 

2 

0 

0 

0 

0 

0 

0 

5.     FLEET  UNIT 

12 

6 

6 

6 

6 

6 

6 

6 

6 

9 

9 

9 

6.     SHORE  DUTY 

0 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

TABLE  3.  CURRENT  INCUMBENTS 

ACTIVITIES/TOURS 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

1.     POSTGR  EDUC 

0 

124 

133 

12 

12 

28 

8 

7 

2 

0 

2 

0 

2.     JPME 

0 

0 

3 

0 

2 

11 

9 

12 

15 

7 

4 

1 

3.     JOINT  TOUR 

0 

1 

3 

3 

3 

12 

15 

36 

91 

41 

39 

19 

4.     SWO  EDUC 

0 

24 

55 

69 

15 

9 

0 

0 

0 

0 

0 

0 

5.     FLEET  UNIT 

3744 

765 

168 

299 

465 

253 

225 

194 

313 

172 

64 

46 

6.     SHORE  DUTY 

0 

466 

530 

126 

88 

119 

117 

244 

428 

252 

137 

83 

TABLE  4.  CURRENT  ACCESSIONS 

ACTIVITIESfTOURS 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

1.     POSTGR  EDUC 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2.     JPME 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3.     JOINT  TOUR 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

4.     SWO  EDUC 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

5.     FLEET  UNIT 

315 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

6.     SHORE  DUTY 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

23 


TABLE  5.  CURRENT  TRANSFER  RATES 


CURRENT  TRANSFER  RATES  WHEN 
LEAVING  TOUR  NUMBER  1 


CURRENT  TRANSFER  RATES  WHEN 
LEAVING  TOUR  NUMBER  7 


ACTIVITIES/ACTIVITIES 

1 

2 

3 

4 

5 

6 

1.          POSTGR  EDUC 

0 

0 

0 

0 

0 

0 

2         JPME 

0 

0 

0 

0 

0 

0 

1         JOINT  TOUR 

0 

0 

0 

0 

0 

0 

4         SWO  EDUC 

0 

0 

0 

0 

0 

0 

5.         FLEET  UNIT 

0.05 

0 

0 

0.05 

040 

0J8 

k.          SHORE  DUTY 

0 

0 

0 

0 

0 

0 

CURRENT  TRANSFER  RATES  WHEN 
LEAVING  TOUR  NUMBER  2 

AcnvrriES/AcrrvrriES 

1 

2 

3 

4 

5 

6 

L         POSTGR  EDUC 

0 

0 

0 

1.00 

0 

0 

2         JPME 

0 

0 

0 

0 

0 

0 

1         JOINT  TOUR 

0 

0 

0 

0 

0 

0 

4         SWO  EDUC 

0 

0 

0 

0 

10) 

0 

S         FLEET  UNIT 

0.13 

0 

0 

0 

0.10 

052 

6.         SHORE  DUTY 

0 

0 

0 

024 

0 

0 

CURRENT  TRANSFER  RATES  WHEN 
LEAVING  TOUR  NUMBER  3 

AcnvmES/AcrrvnTEs 

1 

2 

3 

4 

5 

6 

1.          POSTGR  EDUC 

0 

0 

0 

uoo 

0 

0 

2          JPME 

0 

0 

0 

0 

0 

0 

1          JOINT  TOUR 

0 

0 

0 

0 

0 

0 

4          SWO  EDUC 

0 

0 

0 

0 

1.00 

0 

S         FLEET  UMl 

0 

0 

0 

0 

1.00 

0 

6.          SHORE  DUTY 

0 

0 

0 

0i3 

0 

025 

CURRENT  TRANSFER  RATES  WHEN 
LEAVING  TOUR  NUMBER  4 

AcnvmES/ACTivrriES 

1 

2 

3 

4 

5 

6 

L         POSTGR  EDUC 

0 

0 

0 

t) 

0 

0 

2         JPME 

0 

0 

0 

0 

0 

0 

1         JOINT  TOUR 

0 

0 

0 

0 

0 

0 

4         SWO  EDUC 

0 

0 

0 

0 

1j00 

0 

S         FLEET  UNIT 

0.05 

0 

0 

0 

0.80 

0.15 

&         SHORE  DUTY 

(1 

0 

0 

0.50 

0 

0 

CURRENT  TRANSFER  RATES  WHEN 
LEAVING  TOUR  NUMBER  5 

ACTIVITIES/ACTIVITIES 

1 

2 

3 

4 

3 

t 

1.          POSTGR  EDUC 

0 

0 

0 

0 

1.00 

0 

2          JPME 

0 

0 

0 

0 

0 

0 

i         JOINT  TOUR 

0 

0 

0 

0 

0 

0 

4          SWO  EDUC 

0 

0 

0 

0 

101 

0 

5.         FLEET  UNIT 

0.05 

0.07 

0X32 

0 

oas 

0.4O 

h         SHORE  DUTY 

0 

(! 

0 

0 

0.80 

0 

CURRENT  TRANSFER  RATES  WHEN 
LEAVING  TOUR  NUMBER  6 

AcnvrriES/ACTivrriES 

1 

2 

3 

4 

5 

6 

1.          POSTGR  EDUC 

0 

0 

0 

0 

1.00 

0 

2         JPME 

(1 

0 

040 

0 

0j60 

0 

a         JOINT  TOUR 

0 

0 

0 

0 

0.85 

0 

4          SWO  EDUC 

0 

0 

0 

0 

0 

0 

S          FLEET  UNIT 

0.Q5 

0.10 

0.04 

0 

040 

0  40 

6.          SHORE  DUTY 

0 

0 

0 

0 

0.85 

0 

Acn  vrriES/ACTrv  rriES 

1 

2 

3 

4 

3 

6 

L         POSTGR  EDUC 

0 

0 

0 

0 

1j00 

0 

2         JPME 

0 

0 

0.50 

0 

050 

0 

a         JOINT  TOUR 

0 

0 

0 

0 

UOO 

0 

4         SWO  EDUC 

0 

0 

0 

0 

o 

0 

S         FLEET  UNIT 

0.05 

0.15 

0j05 

0 

0.25 

050 

6.         SHORE  DUTY 

0 

0 

0 

0 

UOO 

0 

CURRENT  TRANSFER  RATES  WHEN 
LEAVING  TOUR  NUMBER  8 

AcnvrriES/AcnvrnEs 

1 

2 

3 

4 

5 

6 

1.          POSTGR  EDUC 

0 

0 

0 

0 

IjOO 

0 

2         JPME 

0 

0 

0.90 

0 

0.10 

0 

a         JOINT  TOUR 

0 

0 

0 

0 

103 

0 

4         SWO  EDUC 

0 

0 

0 

0 

0 

0 

5          FLEET  UNIT 

0 

0.15 

0.15 

0 

025 

0A5 

6          SHORE  DUTY 

0 

0 

0 

0 

0A5 

oas 

CURRENT  TRANSFER  RATES  WHEN 
LEAVING  TOUR  NUMBER  9 

ACnVITIES/ACTIVmES 

1 

2 

3 

4 

S 

6 

L         POSTGR  EDUC 

0 

0 

0 

0 

0 

0 

2         JPME 

0 

0 

050 

0 

050 

0 

a         JOINT  TOUR 

0 

0 

0 

0 

103 

0 

4          SWO  EDUC 

0 

0 

0 

0 

0 

0 

5.         FLEET  UNIT 

0 

0j06 

0.02 

0 

0.15 

0j65 

&          SHORE  DUTY 

0 

0 

0 

0 

.45 

025 

CURRENT  TRANSFER  RATES  WHEN 
LEAVING  TOUR  NUMBER  10 

AcnvmES/AcnvmES 

1 

2 

3 

4 

5 

6 

1.         POSTGR  EDUC 

0 

0 

0 

0 

0 

0 

2         JPME 

0 

0 

050 

0 

050 

0 

a         JOINT  TOUR 

0 

0 

0 

0 

0.95 

0 

4         SWO  EDUC 

0 

0 

0 

0 

0 

o 

S          FLEET  UNIT 

0 

0 

0.10 

0 

0.15 

040 

6.         SHORE  DUTY 

0 

0 

0 

0 

0.20 

oao 

CURRENT  TRANSFER  RATES  WHEN 
LEAVING  TOUR  NUMBER  11 

ACTIVITIES/ACTIVITIES 

1 

2 

3 

4 

3 

6 

1.          POSTGR  EDUC 

0 

0 

0 

0 

0 

0 

2          JPME 

0 

0 

0 

0 

0 

0 

a          JOINT  TOUR 

0 

0 

0 

0 

0.80 

0 

4          SWO  EDUC 

0 

0 

0 

0 

0 

0 

5.          FLEET  UNIT 

0 

0 

020 

0 

0 

oas 

&          SHORE  DUTY 

0 

0 

0 

0 

oao 

035 

24 


TABLE  6.  EXPECTED  NUMBERS  OF  OFFICERS  4  QUARTERS  FROM  NOW 


ACTIVITIES/TOURS 

1 

2 

3 

4 

5 

6 

7 

8 

9 

ID 

11 

12 

1.     POSTGR  EDUC 

0 

126 

135 

8 

14 

26 

10 

10 

0 

0 

0 

0 

2.     JPME 

0 

0 

0 

0 

0 

11 

8 

11 

10 

6 

0 

0 

3.     JOINT  TOUR 

0 

0 

0 

0 

0 

14 

24 

37 

84 

37 

35 

14 

4.     SWO  EDUC 

0 

31 

60 

78 

16 

0 

0 

0 

0 

0 

0 

0 

5.     FLEET  UNIT 

3756 

755 

162 

328 

464 

272 

235 

183 

243 

280 

103 

61 

6.     SHORE  DUTY 

0 

457 

530 

130 

74 

185 

127 

200 

317 

317 

137 

75 

TABLE  7.  EXPECTED  NUMBERS  OF  OFFICERS  10  QUARTERS  FROM  NOW 


ACTIVITIES/TOURS 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

1.     POSTGR  EDUC 

0 

125 

131 

0 

16 

23 

14 

12 

0 

0 

0 

0 

2.     JPME 

0 

0 

0 

0 

0 

11 

9 

12 

9 

4 

0 

0 

3.     JOINT  TOUR 

0 

0 

0 

0 

0 

12 

32 

33 

79 

26 

29 

14 

4.     SWO  EDUC 

0 

31 

58 

77 

17 

0 

0 

0 

0 

0 

0 

0 

5.     FLEET  UNIT 

3774 

749 

169 

342 

495 

276 

267 

192 

193 

373 

150 

80 

6.     SHORE  DUTY 

0 

449 

526 

132 

64 

248 

142 

155 

190 

319 

158 

73 

TABLE  8.  EXPECTED  NUMBERS  OF  OFFICERS  40  QUARTERS  FROM  NOW 

ACTIVITIES/TOURS 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

li 

12 

1.     POSTGR  EDUC 

0 

126 

131 

0 

17 

25 

14 

17 

0 

0 

0 

0 

2.     JPME 

0 

0 

0 

0 

0 

12 

10 

17 

13 

4 

0 

0 

3.     JOINT  TOUR 

0 

0 

0 

0 

0 

13 

34 

42 

113 

31 

34 

20 

4.     SWO  EDUC 

0 

32 

59 

76 

16 

0 

0 

0 

0 

0 

0 

0 

5.     FLEET  UNIT 

3780 

756 

170 

346 

505 

289 

340 

255 

229 

311 

132 

73 

7.     SHORE  DUTY 

0 

454 

524 

131 

69 

269 

153 

227 

232 

243 

166 

84 

25 


LIST  OF  REFERENCES 

[1]  Drescher,  Richard  B.,  'The  Effect  of  the  Goldwater-Nichols  Department 
of  Defense  Reorganization  Act  on  Tactical  Aviation  Pilot  and  NFO 
Career  Paths,"  Master's  Thesis,  Naval  Postgraduate  School,  Monterey, 
CA,  September  1989. 

[2]  Feller,  W.,  An  Introduction  to  Probability  Theory  and  Its  Applications, 
Volume  n,  John  Wiley,  1966,  pp.  189-193  and  374-377. 

[3]  Johnson,  Joseph  L.,  Jr.,  "Analysis  of  the  Joint  Reorganization  Act's 
Impact  on  Personnel  Flow  in  the  Surface  Warfare  Officer  Community," 
Master's  Thesis,  Naval  Postgraduate  School,  Monterey,  CA,  March  1989. 

[4]  Milch,  Paul  R.  (1988).  "An  Analytical  Model  for  Forecasting  Navy  Officer 
Career  Paths,"  NPS  Technical  Report  #NPS55-88-009,  September  1988. 


26 


INITIAL  DISTRIBUTION  LIST 

NO.  OF  COPIES 

Library  (Code  0142)  2 

Naval  Postgraduate  School 
Monterey,  CA    93943-5000 

Defense  Technical  Information  Center  2 

Cameron  Station 
Alexandria,  VA    22314 

Office  of  Research  Administration  (Code  012)  1 

Naval  Postgraduate  School 
Monterey,  CA     93943-5000 

CDR  William  D.  Ferree  1 

Officer  Allocation  and  Distributable  Strength  Projection  Branch 

Naval  Military  Personnel  Command 

NMPC-454 

Washington,  DC    20350-2000 

Professor  Peter  Purdue,  Code  OR/Pd  1 

Naval  Postgraduate  School 
Monterey,  CA    93943-5000 

Professor  Paul  R.  Milch,  Code  OR/Mh  20 

Naval  Postgraduate  School 
Monterey,  CA    93943-5000 

Center  for  Naval  Analyses  1 

4401  Ford  Avenue 
Alexandria,  VA    22302-0268 

Library  (Code  OR)  1 

Naval  Postgraduate  School 
Monterey,  CA     93943-5000 

Institute  for  Defense  Analysis  1 

1800  North  Beauregard 
Alexandria,  VA    22311 


27 


Arthur  P.  Hunter,  Jr. 
Professor  and  Chairman 
Dept.  of  Industrial  Engineering 

and  Management  Sciences 
Northwestern  University 
Evanston,  IL    60201-9990 

Koh  Peng  Kong 
OA  Branch,  DSO 
Ministry  of  Defense 
Blk  29  Middlesex  Road 
SINGAPORE  1024 

Commanding  Officer 

ATTN:   Mr.  Murray  Rowe 

Navy  Personnel  Research  and  Development  Center 

San  Diego,  CA    92152 

Commanding  Officer 

ATTN:    Barry  Siegel 

Navy  Personnel  Research  and  Development  Center 

San  Diego,  CA    92152 

Mr.  R.  Morton 

%DSAI 

350  Fortune  Terrace 

Rockville,  MD    20854-2995 

Operations  Research  Center,  Rm.  E40-164 
Massachusetts  Institute  of  Technology 
ATTN:   R.  C.  Larson  and  J.  F.  Shapiro 
Cambridge,  MD    02139 

James  W.  Mosteller,  Code  111 

Naval  Personnel  Research  and  Development  Center 

San  Diego,  CA    92152-6800 

Dr.  M.  J.  M.  Posner,  Professor 
Department  of  Industrial  Engineering 
Faculty  of  Applied  Science  and  Engineering 
University  of  Toronto 
Toronto,  CANADA     M5S  1A4 


28 


LT  Joseph  Johnson 
442  Warely  St. 
Portsmouth,  RI     02871 

Ms.  Annette  Stephens,  2d  Lt,  USAF 

Scientific  Analyst 

Cmd  Analysis  Division 

DCS /Plans  &  Requirements 

Department  of  the  Air  Force 

Headquarters  Air  Training  Command  (ATC) 

Randolph  Air  Force  Base,  TX    78150-5001 


29 


DUDLEY  KNOX  LIBRARY 


mill 

ii 

hi 

;mi 
111, 

] 

biii 

3  2768  00342571  1 


